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Abstract 


This study investigates the variations in body size, bioacoustic traits, 
and cuticular hydrocarbon (CHC) profiles in different populations of the 
bush cricket species Isophya autumnalis Karabag, 1962. Within-population 
body size variations, particularly those associated with distinct habitat dif- 
ferences and climate shifts within their local distribution ranges, suggest 
that ecological factors affect morphological characteristics. Sexual size di- 
morphism (SSD) in I. autumnalis may affect reproductive behavior and 
strategies, potentially influenced by the bioacoustic environment. While 
male calling songs exhibit temporal variations across populations, sug- 
gesting differences among allopatric populations, CHC profiles, known to 
undergo selection under various climatic conditions, also vary noticeably 
across local populations. These findings highlight the importance of un- 
derstanding within-species variations for the conservation of Isophya and 
similar taxa in the face of habitat threats. Overall, this study contributes to 
a comprehensive understanding of how morphology and bioacoustic be- 
havioral traits are shaped over short distances in allopatric populations of 
species with limited mobility, such as I. autumnalis, providing insights into 
adaptation processes and highlighting the urgency of conservation efforts 
for endemic species in Anatolia. 
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Introduction 


There are four main speciation types: allopatric, peripatric, 
parapatric, and sympatric (Howard and Berlocher 1998, Tregenza 
et al. 2000, Butlin et al. 2008, Fitzpatrick et al. 2009). The most 
critical factor accelerating allopatric speciation is topographic 
isolation that prevents reproductive unity between populations. 
Many studies have shown genetic and morphological differentia- 
tion from the ancestral population due to genetic drift and diver- 
gent selection caused by various ecological conditions due to the 


reduction or complete cessation of gene exchange between geo- 
graphically separated local populations of the common ancestral 
population (Hemp et al. 2010, Langerhans and Riesch 2013, Tay- 
lan and Sirin 2016, Sevgili et al. 2018). However, many questions 
remain unanswered. For example, questions regarding the extent 
to which ecological factors (such as climate, vegetation, predator- 
prey relationships, etc.) lead to reproductive isolation (habitat, 
temporal, or gametic) in local populations within a species’ range 
(Sobel et al. 2010, Langerhans and Riesch 2013); the degree of vari- 
ation that may arise in any trait within local populations as they 
adapt to varying ecological conditions; and which traits (mor- 
phology, behavior, etc.) are more plastic if allopatric speciation 
occurs, remain insufficiently addressed. 

Understanding the degree of variation in biochemical, behav- 
ioral, and morphological characters within populations is crucial. 
Insect mating recognition systems often rely on species-specific 
chemical compound profiles driven by strong sexual selection 
(Chapman et al. 2000, Kather and Martin 2012, Moore et al. 2021). 
In contrast, the evolutionary rates of the taxonomic characters 
used to define species are not uniform. Environmental conditions, 
population density, selection pressure, behavioral differences, and 
mutation rates can vary within local populations, influencing the 
speed of change (Bromham 2009, Adams 2013). For example, be- 
havioral characters (social behaviors or bioacoustic characters) can 
evolve more rapidly than morphological traits (Losos 1990, Heller 
2006, Pitchers et al. 2014). Orthoptera is one of the insect groups 
that provides the best evidence to support this claim. 

In particular, research on meadow grasshoppers of the genus 
Chorthippus Fieber, 1852 (Acrididae) shows that whereas its species 
are very difficult to distinguish morphologically, genetically, and 
ecologically, they can be clearly differentiated through species-spe- 
cific bioacoustic traits (Sirin et al. 2010, Nolen et al. 2020, Tarasova 
et al. 2021). In contrast, studies of the bush cricket species of the 
genus Poecilimon Fischer, 1853 have identified significant interspe- 
cies differences in their morphological, genetic, ecological, and be- 
havioral characters (Heller 1988, Sevgili et al. 2008, Kocinski et al. 
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2022). In some bush cricket genera, such as Eupholidoptera Maran, 
1953 and Parapholidoptera Maran, 1953, bioacoustic characters are 
quite similar across species, while morphological and molecular 
characters are more likely to be selected for (Heller 2006, Ciplak 
et al. 2021). Specifically, it is important to study local populations 
of species with short wings and limited mobility to determine in- 
traspecific variation in characters because these local populations 
have a fragmented distribution, so there are only narrow zones 
where gene exchange can occur. In populations with limited mo- 
bility due to geographical isolation, selection for pre-breeding iso- 
lation mechanisms can be strong, thereby accelerating allopatric 
and parapatric speciation (Eweleit et al. 2015, Sevgili 2018). 

Turkey has a highly diverse geography and rich biological di- 
versity with different genetic patterns across short distances. In the 
present study, we investigated the variability of several taxonomic 
characters in local subpopulations of Isophya autumnalis Karabag, 
1962 (JA). The genus Isophya Brunner von Wattenwyl, 1878 is dis- 
tributed in the Palaearctic region but is especially concentrated in 
the Balkans and Anatolia and includes species with short wings 
and limited mobility. Out of 90 species globally, the 44 species 
and subspecies distributed in Anatolia include 77% local endemic 
species with narrow distributions (Sevgili 2004, 2018, 2020, Sev- 
gili et al. 2006, 2012). We tried to reach other potential local IA 
populations previously reported from the type locality, namely the 
Zigana Mountains, Trabzon. We hypothesized that there are dif- 
ferences in morphological/morphometric, bioacoustic, and bio- 
chemical (cuticular hydrocarbons, CHC) characters between IA 
subpopulations at localities that are somewhat distant and have 
different climatic conditions and vegetation structures. Body size, 
which is a vital trait for organisms, is highly variable among local 
populations because of its high sensitivity to changes in biotic and 
abiotic factors (Chown and Gaston 2010). In addition to having 
different distribution patterns, local population body sizes may 
also vary due to different climate (Bidau and Marti 2007, Eweleit 
and Reinhold 2014), food/water resources (Stillwell et al. 2007), 
habitat type/width and vegetation structure (K6nig and Krauss 
2019), size-dependent predation pressure (Remmel and Tammaru 
2009), population density (White et al. 2007), and general life his- 
tory traits (Roff 2002, Chown and Gaston 2010). 

Isophya’s distinctive species-level taxonomic characters derive 
from morphology, bioacoustics (Heller et al. 2004, Sevgili 2004, 
Sevgili et al. 2006, Unal 2010), mitochondrial (cyt b, COI, COII, 
and ND2), and nuclear (ITS1 and ITS2) molecular markers 
(Grzywacz-Gibala et al. 2010, Chobanov et al. 2017). However, 
the molecular findings and taxonomy based on morphology and 
bioacoustic data (species-specific male calls) do not fully align 
regarding the phylogeny and differentiation of Isophya species 
(Grzywacz-Gibala et al. 2010, Chobanov et al. 2017). JA belongs 
to the I. zernovi species group (Sevgili 2004, 2020, unpublished 
data Sevgili et al.). While ecological differentiation plays a crucial 
role in speciation, it is important to recognize that sexual selection 
can also interact with it, further influencing the process (Weissing 
et al. 2011). 

CHCs, which consist of an intricate blend of unsaturated hy- 
drocarbons, methyl-branched alkanes, and n-alkanes secreted from 
arthropod cuticles, are a biochemical taxonomic character in social 
insects (Martin et al. 2008, Kather and Martin 2012, Pena-Carrillo 
et al. 2021) and some other insect groups (Chorthippus sp., Finck 
et al. 2016; Drosophila sp., Davis et al. 2021; ladybugs, Redjdal et 
al. 2023). CHCs have many functions. In terrestrial environments, 
they primarily prevent dehydration while also serving as sexual 
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signal molecules in a wide variety of chemical communication sys- 
tems (Drijfhout et al. 2010, Ingleby 2015, Holze et al. 2021). 

Rather than addressing these functions, we investigated the 
taxonomic utility of CHC profiles based on their variation across 
local populations of JA. However, we also hypothesize that CHC 
profiles in Isophya populations will be less variable than morpho- 
logical/morphometric and bioacoustic characters. This is because 
Isophya species are not eusocial insects, and male calling songs are 
used for communication, such as finding mates. We also hypoth- 
esized that males and females within a population would have 
dimorphic CHC profiles. The results of this study are important 
for understanding which characters are particularly important for 
speciation processes in narrowly distributed and endemic popu- 
lations. The results can also aid efforts to conserve narrowly dis- 
tributed endemic species by clarifying which traits are more vari- 
able in allopatric populations. Given that Isophya species have a 
fragmented distribution pattern due to their limited mobility and 
unique habitat preferences (Nuhli¢kova et al. 2023), and that 
gene exchange is largely restricted, it is important to determine 
how much the multiple taxonomical characters differ between lo- 
cal populations. 


Materials and methods 


The plump bush crickets. —The model for our study, JA, is an en- 
demic species described by Karaba& in 1962 based on mate- 
rial collected during a field study in September from the Zigana 
Mountains (Trabzon, Turkey). In species of the genus Isophya, 
the wings of both sexes are shortened and do not reach half of 
the abdomen, and the wings, especially in males, have a species- 
specific sound-producing organ (stridulatory organ) (Heller and 
von Helversen 1986, Sevgili et al. 2012, Iorgu et al. 2013). All 
Isophya species are plump-green bush-crickets with limited mo- 
bility (Sevgili 2004, Nuhliékova et al. 2023). They are generally 
among the earliest bush-cricket species to emerge in the field at 
the end of winter. They hatch earlier (around February) in the 
southern parts of the genus’s distribution range (e.g., Southeast- 
ern Anatolia, Isophya sikorai Ramme, 1951 and later (around June) 
in northern Anatolia and the high mountains (e.g., Mountains 
in Black Sea Region, I. rizeensis Sevgili, 2003 (Sevgili 2004). Field 
observations show that males hatch earlier than females. Isophya 
species generally reach adulthood after five molts during the nym- 
phal stage (Can 1958). Males first produce a calling song at 2-3 
days old, while successful mating generally begins after the sev- 
enth day (when sperm transfer occurs) (Uma and Sevgili 2015). 
IA is a typical mountain species, and the northern population is 
distributed within alpine zones at high altitudes in Zigana Moun- 
tain. Nymphs can be found in suitable habitats in June and adults 
until mid-September. The species was known only from the type 
locality, but we found populations further south as a result of our 
field studies. We observed several morphological and phenologi- 
cal differences in local populations and investigated the variation 
in taxonomic characters. 


Spatial differences of the sampling locations.—Species’ presence in a 
habitat, their ecological niche, and distribution patterns, includ- 
ing habitat quality, seasonality, and life cycles, can often be un- 
derstood by utilizing remote sensing variables (Leitao and Santos 
2019, Schwager and Berg 2021). In this study, we determined some 
environmental characteristics of the five sampled localities using 
variables provided by satellite remote sensing (S-RS) products, 
which are used to measure species—environment relationships 


JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) 


E. KIRAN OZDEMIR, H. SEVGILI AND E. BAGDATLI 


that predict species distributions. In particular, it was found that 
S-RS products provide more reliable spatial resolution when used 
together with bioclimatic data to distinguish the habitat charac- 
teristics of species (Pinto-Ledezma and Cavender-Bares 2020). 
Normalized difference vegetation index (NDVI), which measures 
green vegetation cover, and normalized difference moisture index 
(NDMI), which monitors vegetation water content and drought, 
were used to compare the locations studied. NDVI, which ranges 
from -1 to 1, isa measure of the health of vegetation based on how 
plants reflect light at certain wavelengths. Values close to zero in- 
dicate barren areas, such as rocky and sandy areas, while a positive 
1 indicates forested areas with mild temperatures and abundant 
rainfall. NDMI also ranges from -1 to 1 and provides important 
data on the water content and drought status of vegetation. Nega- 
tive values of NDMI correspond to barren soil, values around zero 
indicate water stress, and high positive values represent canopy 
without water stress. 

For this purpose, the Sentinel-2B level-2A product was down- 
loaded using the European Space Agency (ESA) data portal 
(https://browser.dataspace.copernicus.eu/). Satellite images were 
downloaded after applying a search constraint to ensure cloud 
cover was below 5%. The files were downloaded with high resolu- 
tion selected. All viewing and analysis processes were conducted 
using the “raster” and “ggplot2” R packages. Satellite data from 
July 15, 2018 to August 31, 2018 were evaluated, which is the pe- 
riod when the adult stages of IA passed through and samples were 
collected. NDVI and NDMI data were downloaded for both July 
and August for a snow-free and cloud-free 2 km? area over the loca- 
tion where the specimens were collected. In addition, Sentinel -2 
bands BO4 (red) and BO8 (NIR-near-infrared) were used to calcu- 
late the NDVI index. The normalized difference vegetation index, 
NDVI = (NIR-RED)/(NIR+RED) equation was used to calculate 
the normalized difference moisture index, NDMI = (NIR-B11)/ 
(NIR+B11) (Piedelobo et al. 2019, Lastovicka et al. 2020). B11 and 
12 short-wave infrared (SWIR) bands are often used for examining 
surface features such as vegetation, soil, and water. Calculations 
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and plots were performed in the R environment using raster and 
rasterVis packages (Hijmans and van Etten 2012, Perpifan and Hi- 
jmans 2023). SWIR band measurements, a technique for classify- 
ing species, utilizes either the continuous spectrum from visible 
to shortwave infrared, specific spectral curve shape-based features, 
or vegetation indices. Linking spectral data with biochemical and 
physiological ancillary data can enhance classification accuracy and 
relate species’ spectral separability to their characteristic ecophysi- 
ological traits. In this context, vegetation indices and spectral fea- 
tures can be considered semi-quantitative biochemical parameters 
of the reflectance spectrum, aiding in assessing the physiological 
status of vegetation (Grofe-Stoltenberg et al. 2016). Monthly/aver- 
age annual temperature and relative humidity data for the studied 
regions by single point method were obtained from Prediction of 
Worldwide Energy Resources (POWER) (https://power.larc.nasa. 
gov/beta/data-access-viewer/) by filtering for “temperature at 2 
meters” and “relative humidity at 2 meters.” A table was created by 
adding one year before and after the sampling year (2018). Since 
no specific data for Pekitn Mountain could be produced from this 
website, it was not included in the table. 


Field study.—During field studies conducted in Trabzon and 
Giimushane provinces in June-July 2018, populations were found 
at five different locations in, approximately, a line from north to 
south (Fig. 1). From each location, 60 male and female nymphs of 
the last stage were collected and transported live to the laboratory 
in 30 x 20 x 20 cm cages with a piece of plant from the habitats 
where they were collected for feeding (Suppl. material 1: fig. $1). 


Maintenance of the bushcrickets.—In the laboratory, the samples 
were transferred to larger 50 x 30 x 30 cm cages, with males and 
females separated. The cages were cleaned every two days and sup- 
plied with the same amount of nettle (Urtica dioica), Taraxacum 
sp., lettuce, and apple slices. The adult males and females were 
transferred to smaller individual cages and given the above-men- 
tioned food. Water was sprayed into the cages every day to provide 
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Fig. 1. Map of sampling localities of Isophya autumnalis. The sampling sites are also the latest known distribution points of the species. 
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moisture. Adult males began producing sounds on the second or 
third day, and calling songs were recorded after the third day. Spec- 
imens were collected from field locations at different times due to 
differences in aspect, vegetation structure, and distance between 
locations. However, feeding conditions were standardized, and 
care was taken to ensure that the specimens were of the same age. 


Bioacoustics.—To record the males’ two-syllable calling songs, 
each male was placed in a small cage and taken to an isolated 
room where it could not hear the calling songs of other males. 
Songs were recorded from each individual for at least 3 minutes. 
All recordings were made under the same temperature conditions 
(25-27°C). For song recordings, a condenser ultrasonic micro- 
phone (Avisoft Bioacoustics CM16/CMPA, sampling rate 96 kHz) 
connected to a digital recorder, TASCAM HD-P2, was used. The 
recorded songs were analyzed using Raven Pro software (Raven 
Sound Analysis 2022) by taking the necessary spectral measure- 
ments. The males were then preserved at -80°C together with the 
females (all 8 days old) based on their field locations to examine 
the CHC profiles and morphological measurements. Song meth- 
odology and terminology follows Heller et al. (2004), Sevgili et al. 
(2006) and Sevgili et al. (2012). 


Morphology/morphometry.—The taxonomically important anatomi- 
cal parts of each specimen, such as the head, fastigium, pronotum, 
wing, hind femur, cercus, subgenital plate, epiproct, and oviposi- 
tor of females, were photographed with a stereo zoom microscope 
(Suppl. material 1: fig. $2). The morphology, number of teeth, and 
shape of the male stridulatory organs differ significantly between 
species (Sevgili 2004). Therefore, to identify variations in the 
stridulatory organs between local populations within the species, 
photographs of the organ on the upper left wing of the males from 
each population were taken using a scanning electron microscope 
(SEM Hitachi, SU-1510). Morphometric measurements of all pho- 
tographs were taken using the ImageJ program (https://image). 
nih.gov/ij/download.html). In measurements of the stridulatory 
file, the focus was on the shortest distance of the file, the distance 
to the wing proximal and distal edges, and the number of teeth 
(Suppl. material 1: fig. S3). 

Regarding morphology, measurements were taken of the fas- 
tigium length, fastigium base width, interocular distance, dorsal 
pronotum length, pronotum width, lateral pronotum length, right 
posterior femur length, epiproct length and width, subgenital 
plate length, shortest proximal and distal cercus lengths in males, 
and length of the ovipositor and subgenital plate of females (see 
Suppl. material 1: fig. $2). 

To identify potential size dimorphism in local populations, 
sexual size dimorphism (SSD) was estimated for each population 
based on the following formula: SSD = (size of the larger sex/size 
of the smaller sex). A positive result indicates that females are the 
larger sex and vice versa (Lovich and Gibbons 1992, Stillwell et 
al. 2007). The hind femur length was used in the formula because 
it is the best predictor of body size in bush-crickets (Eweleit and 
Reinhold 2014). 


CHC extraction and analysis.—Determination of the specimens’ 
CHC profiles followed the methods traditionally used for research 
on Orthoptera (Neems and Butlin 1994, Thomas and Simmons 
2008, 2010, Thomas et al. 2011, Veltsos et al. 2012, Steiger et al. 
2013). CHC data were obtained using a Shimadzu GCMS QP2010 
Ultra gas chromatograph-mass spectrometer (GC/MS). Individu- 
als previously frozen at -80°C were kept in a desiccator for 45 
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minutes to reach room temperature and to get rid of excess wa- 
ter. Each individual was then placed in its entirety in a separate 
conical-based glass centrifuge tube and 4 ml, 10 ppm pentadecane 
(internal standard) containing hexane was added. After waiting 
for 30 seconds in hexane and 5 minutes including vortex, the up- 
per layer of the solution was taken from the vial using dispos- 
able glass pipettes and placed in the GC/MS device for analysis. 
To determine the CHC profiles in the bush-crickets, each sample 
(1 pl) was injected into the device for split-mode (10:1) analysis 
in a Stabilwax column (Restek/Rtx-5) with an inner diameter of 
30 m x 0.25 mm and chromatographic purity under helium gas. 
The analysis started at 50°C for 1 min, and then the temperature 
was raised to 250°C at 15°C per min. The final temperature was 
increased to 320°C at 3°C per min, and the sample was held at 
this temperature for 5 min. The extract was separated using the 
temperature difference in the column. The transfer of compo- 
nents to GC/MS was done at 250°C. Hexane was added to the 
analysis every day to prevent contamination. Peak numbers were 
coded taking into account the retention time for data analysis. The 
names of basic hydrocarbons were identified by examining GC/ 
MS libraries (NIST-17 and Willey libraries) and other published 
data. The fragmentation profiles were acquired and analyzed with 
theMS and MSD ChemStation software. 


Statistics. —Principal component analysis (PCA) was applied to the 
measurement data from the morphological parts for both sexes 
(66 females; 57 males), as shown in Suppl. material 1: fig. S2. 
Cercus and ovipositor measurements were included for males and 
females, respectively. 

PCA was applied to number of teeth, distal distance, proximal 
distance, and stridulatory file length, taken from the SEM image of 
the stridulatory organs of the male wings (Fig. 6; Supp]. material 
1: fig. $3). All components with eigenvalues greater than 1 were 
included in the statistical analysis. The CHC phenotypes of the 66 
females and 57 males collected from all localities were compared. 
Hydrocarbons were identified based on retention times for their 
mass spectra. The CHC data were evaluated statistically follow- 
ing the methods of related studies (Thomas and Simmons 2009, 
Steiger et al. 2013, Finck et al. 2016). 

Before the analysis, the area under each peak on the chroma- 
tograph was divided by the area of the internal standard (penta- 
decane) to control for variations in CHC extraction efficiency be- 
tween specimens. These ratios were then transformed using log10 
to ensure that the data were normally distributed. The data for 
each peak were then analyzed by PCA using the “Factoextra” and 
“FactoMineR” R packages (Lé et al. 2008, Kassambra and Mundt 
2020). The first two components (PC1, PC2) were compared 
based on locality and sex and visualized using the “ggplot2” R 
package (Wickham 2016). The first five PC scores with eigenvalues 
greater than 1 (PC1) were analyzed multivariately using the de- 
pendent variables of locality and sex and the interaction variable 
sex* localities. Tukey’s HSD post hoc analysis was applied to com- 
pare each PC score between localities. All analyses were performed 
in R Studio (Rstudio Team 2021). 


Results 


Spatial differences of the localities. —The Zigana and Krom Valley re- 
gions were evaluated together due to their geographical proximity 
and the availability of similar climate data. All the regions evalu- 
ated exhibited different patterns of temperature, humidity, and 
precipitation (no data for Pekiin) (Suppl. material 2: table $1). 
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The pixel values of the July and August plots of NDVI, NDMI, 
and SWIR data vary across localities (Suppl. material 1: figs S4-S8; 
Suppl. material 2: table $2). The mean NDVI, NDMI, and SWIR 
data for the locations are also different in both months. Zigana 
has very low mean values for all three data (for both months), 
while Vauk and Péske have the highest values (Suppl. material 1: 
figs S9-S11). 


Morphometry.—Suppl. material 2: table S3 provides the mean and 
standard deviation values for various morphological and bio- 
acoustic measurements of the species. Regarding body size (spe- 
cifically hind femur length), males from the Zigana and Krom 
Valley populations were significantly larger than those from the 
southern populations (ANOVA, F,., = 23.72, p < 0.001). How- 
ever, there were no significant body size differences between the 
other groups. For females, the pattern was different. Only the Zi- 
gana population exhibited larger body size (ANOVA; F, ,, = 12, 
p < 0.001). There were no significant body size differences among 
the other female populations. 

Regarding the PCA of the morphometric measurements of 
males, the first four components, which collectively explained 
67.67% of the variance, all had eigenvalues greater than 1. 
Comparing populations based on PC1, the Zigana population 
was significantly distinct from the other populations (ANOVA, 
F, ,, = 21.15, p< 0.001) (Fig. 3). There were no significant differenc- 
es between the mountain populations of Krom, Vauk, and Péske 
(Dien ae opi = 01226, Dice... = 0.363, Ph poe = 10.066), whereas 
the Vauk and Pekiin populations were distinct (Py, ,. pexan = 9-029). 

PC1 was predicted by the following variables in descending 
order of importance: interocular distance (15.72%), hind femur 
length (14.46%), epiproct width (14.46%), lateral pronotum 
length (13.77%), and cercus length (13.04%). The PC2 compo- 
nent was significantly influenced by fastigium length (30.16%), 
proximal width of fastigium (38.11%), and epiproct length 
(26.86%). 

Regarding the PCA of the female morphological measure- 
ments, the first three components, collectively explaining 68.72% 
of the variance, had eigenvalues greater than 1. Comparing the 
populations based on PC1 and PC2, only the Zigana population 
was significantly distinct from the others (ANOVA; F, ,, = 28.91, 
p < 0.001) (Fig. 4). 

In terms of common morphological structures measured in 
both males and females (fastigium, interocular distance, prono- 
tum, hind femur, and subgenital plate), PC1, which explained 
45.49% of the variance, exhibited clear sexual dimorphism, with 
females being significantly larger (t-test; t = -17.71, df = 119.91, 
p < 0.001). However, for pronotum length, which is an impor- 
tant indicator of body size, like hind femur length in bush crickets 
(Eweleit and Reinhold 2014), there was no significant size differ- 
ence between males and females (t-test, t = 0.924, df = 99.492, 
p = 0.358). 


Sexual size dimorphism.—Although there was a decreasing SSD 
trend from north to south in the local IA populations, this dif- 
ference was not statistically significant (ANOVA; F, ., = 1.897, 
t = -1.377, R? = 0.034, p = 0.174) (Fig. 5). 


Male calling song.—The main syllable typically exhibits a crescen- 
do structure, transitioning to a decrescendo slightly past the mid- 
point. Sometimes, an “after-click” impulse is created at the end of 
the syllables (Fig. 2). Temporal measurements of the songs were 
made as shown in Fig. 2. 
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Stridulatory organ and bioacoustics.—SEM photos of stridulatory 
files of the populations are given in Fig. 6. PC1 was most deter- 
mined by stridulatory organ length (46.36%), followed by proxi- 
mal distance (27.92%), number of teeth in the stridulatory organ 
(16.16%), and distal distance (9.55%). PC2 was most determined 
by number of teeth (37.58%), followed by right proximal distance 
(36.36%), left proximal distance (19.55%), and shortest distance 
length of stridulatory organ (6.50%). 

Stridulatory organ length was significantly positively corre- 
lated with number of teeth (R* ,,, = 0.33, p < 0.001). Regarding 
the number of teeth, the Krom Valley population had signifi- 
cantly fewer than the others except the Péske Dagi population 
(ANOVA, F, ,, = 8.419, p < 0.001). There were no significant dif- 
ferences between the other populations. The average number of 
teeth across all populations was 156+13.11 (meanssd) (range: 
128 to 185). 

Regarding the PCA analysis of the male stridulatory organ 
measurements, the first two components, which explained 80% 
of the variance, both had eigenvalues greater than 1. Com- 
paring populations based on PC1, there were significant dif- 
ferences (ANOVA; F,,, = 13.47, p < 0.001) (Fig. 7). Post-hoc 
analysis (TukeyHSD) revealed that, similar to the other mor- 
phological measurements, the Zigana population differed 
significantly from all the others (P,,...-vauk = 9-008, Pricana-poske 
tOO0U Ds, caper = OCOD yi a troars 0.001). Among the oth- 
ers, only the Péske and Pekiin populations differed significantly 
(p = 0.002). 

PCA was also conducted for 12 variables measuring the tempo- 
ral features of the male calling songs. These included five temporal 
variables and seven impulse-related variables during crescendos 
and decrescendos, and second syllable impulse counts. The first 
four PCs, explaining 77.65% of the variance, all had eigenvalues 
greater than 1. PC1, which explained the most variance (32.29%), 
differed significantly across populations (ANOVA, F, ,,,. = 919.9, 
p < 0.001). Post-hoc analysis indicated significant differences be- 
tween all populations except the Vauk and Péske Mountain popu- 
lations (p = 0.396) (Fig. 8). 

The variables contributing the most to PC1 were number of 
impulses and crescendos in the first syllable (13.04% and 17.64%, 
respectively), the duration of first syllable (11.63%), and the dura- 
tion between the two syllables (15.52%). The variables contrib- 
uting most to PC2 were total duration of the sound (23.66%), 
duration of the first syllable (11.06%), duration between the two 
syllables (8.56%), duration of song period (12.20%), and number 
of impulses in the first syllable (13.04%). 


Cuticular hydrocarbon analysis.—The CHC profile analysis iden- 
tified four types of hydrocarbon structures: alkanes, monome- 
thyl alkanes, dimethyl alkanes, and alkadienes. The following 
straight-chain alkanes were the most frequently represented: 
hexatriacontane (C36 — 13.55%), tritriacontane (C33 - 8.18%), 
hentriacontane (C31 - 7.10%), nonacosane (C29 - 6.61%), tetra- 
triacontane (C34 - 6.39%), heptacosane (C27 - 6.61%), and 
octacosane (C28 - 5.20%). Of the other peaks, the major com- 
ponents were 3,11-dimethylnonacosane (7.21%) and 2-methyl- 
hexacosane (7.10%). 

GC/MS showed that the CHC profile of IA contained differ- 
ent number of compounds for each sex (Figs 9, 10). Except for 
the Zigana population (16 each), females and males in each lo- 
cality had different numbers of CHC components: 18 and 19 in 
Krom Valley; 11 and 17 in Pekiin Mountain; 14 and 19 in Poéske 
Mountain; 11 and 22 in Vauk Mountain. Certain alkanes were 
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Fig. 2. Schematic oscillograms showing the studied male calling song characters and measurement points. 
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significant differences. Those with the same letters do not have a statistically significant difference from each other. 


found exclusively in each sex: docosane (C22) and dotriacon- 
tane (C32) in females; 3,11-dimethyloctacosane and unidenti- 
fied 26, 27, 28, and 37 carbon atom peaks in males. Further- 
more, small proportions of alkene structures were detected in 
males (1%) but not females. The proportions of alkane struc- 
tures were 71% in males and 82% in females. The proportion 
of branched monomethyl alkane profiles was higher in males 
(17%) than females (5%), whereas the proportion of branched 
dimethyl alkanes was the same (7%). Alkadiene structures were 
barely detected in either sex (1%). 

PC1, PC2, and PC3 clearly separated the local populations, 
while PC1 and PC2 separated individuals according to sex and 
interaction of sex and locality (Table 1, Fig. 11). However, this 
difference was observed only between the of Vauk Mountain 
population and other local populations. Significant differences 
were found in terms of the PC3 component between other loca- 
tions (Table 1). 
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Fig. 5. Sexual size dimorphism in Isophya autumnalis based on 
hind femur length. 
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Table 1. Statistics of the cuticular hydrocarbon variation for five local populations of Isophya autumnalis. 


Effect PCl1 PC2 PC3 PC4 PC5 
ES 113 P ES 113 F, 113 P E 113 P re 113 P 
Locality 8.703 <0.001 pa Id <0.001 4.702 0.001 0.619 0.65 07592 0.669 
ae Ps Peas ee ae 
Sex 162335 <0.001 daeod oA <0.001 0299 0.585 0.486 0.487 0.942 0.333 
Sex:Locality 11.886 <0.001 10.133 <0.001 1.301 0.274 2522 0.045 2.204 0.073 
Tukey's HSD post-hoc test os on a De a 
Poéske M.-Vauk M. <0.004 0.011 O77 0.996 0.993 
Zigana-Vauk M. <0.001 0.204 0.011 0:995 0:995 
Krom Valley-Vauk M. <0.001 07329 0.997 07795 0,932 
Pekiin M.-Vauk M. <0.001 <0.001 0.993 0.981 O21 
Zigana-Péske M. 0.932 0.887 0.002 0.949 1 
Krom Valley-Péske M. 0.741 0-732 0.907 0.586 0.762 
Pekiin M.-Péske M. 0.642 0.892 0.999 0.883 0.834 
Krom Valley-Zigana 0.995 0.999 0.043 0.96 0.805 
Pekiin M.-Zigana 0.99 0.383 0.002 0:999 0.871 
Krom Valley-Pektin M. 0099 0.218 0.949 0.971 0.999 
Discussion 


Péske M. 
Fig. 6. Scanning electron microscope (SEM) images showing ex- 
amples of the male stridulatory file organ by location. 


Body size.—As in many animal groups, body size (measurements/ 
mass) is one of the most important traits directly or indirectly 
affecting physiological processes and adaptability in Orthoptera 
species (Roff 2002, Whitman 2008). Many studies have inves- 
tigated the variation of body size in populations living in areas 
where abiotic and biotic factors are influenced and consequently 
changed by different environmental conditions (Whitman 2008, 
Eweleit and Reinhold 2014). The fragmentation of a species’ dis- 
tribution range exposes local populations to varying environmen- 
tal conditions. This variation can lead to changes in physiologi- 
cal constraints, with increased or decreased spatial and temporal 
habitat heterogeneity playing a potential role. For example, studies 
have shown a decrease in bee body size variation (heterogeneity) 
at higher altitudes (Classen et al. 2017). However, the interplay 
between these environmental and physiological factors is complex 
and context-dependent. The physiological constraints in high-alti- 
tude and harsh environments may filter out extreme phenotypes. 
Thus, the lack of significant variability in body size variables in IA 
populations other than the Zigana population may indicate a lack 
of directional selection in local populations exposed to similar al- 
titudes and ecological conditions. Because JA cannot fly or jump 
long distances, they are expected to be more sensitive to specific 
environmental conditions and exhibit morphological, genetic, 
and behavioral diversity (i.e, Amara alpina (Paykull, 1790) (Car- 
abidae ground beetle), Beckers et al. 2020). Macropterous morphs 
disperse over longer distances than brachypterous morphs (Ponia- 
towski and Fartmann 2011, Nuhlitkova et al. 2023). This model is 
supported by the weak body size variability within local popula- 
tions of JA, with the exception of the population in Zigana Moun- 
tain, which is found in habitats with north-facing slopes and typi- 
cal Black Sea high mountain climatic characteristics. 

An analysis of over 1,500 Orthoptera species indicated a gen- 
eral female-biased body size dimorphism, which is stronger in 
Caelifera than in Ensifera (Hochkirch and Gréning 2008). The 
PC1 obtained from the measured body parts of all our IA popula- 
tions revealed a pattern consistent with this analysis, although no 
significant sex differences were found in hind femur lengths. This 
suggests that SSD, in terms of general body size parameters, can- 
not always be determined based on a single character. The larger 
body size of females is usually an evolved response to intraspecific 
competition focused on reproductive success whereby increased 


JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) 


E. KIRAN OZDEMIR, H. SEVGILI AND E. BAGDATLI 


A Stridulatory Organ 


4 


PC2(39.16% 


1 
be 


241 


Localities 


— 


PCI (41.83%) 


1 
be 


Vauk M. 


Zigana M. M 
Localities 


Krom Valley 


Pekiin M. 


Poske M. 


Fig. 7. Principal component analysis (PCA) plot based on the variables proximal-distal shortest length, distance to the distal and 
proximal edges, and number of sound teeth of male stridulatory organs of Isophya autumnalis (A). Boxplot of the measurements among 
different populations (B). Different letters on top of boxplots indicate significant differences. Those with the same letters do not have 


a Statistically significant difference from each other. 


body size is associated with more eggs and reproductive success 
(Hochkirch and Gréning 2008). As with other Isophya and Poecili- 
mon species, IA males are more mobile and therefore spend more 
energy searching for females who respond bioacoustically in order 
to mate (McCartney et al. 2008). SSD and intraspecific body size 
variation can also determine male reproductive behavior strategies, 
thereby shaping certain immune system parameters, such as the 
transfer of more sperm and larger spermatophores to larger females 
(Uma and Sevgili 2015, Sevgili 2022). In the wing-dimorphic bush- 
cricket Metrioptera roeselii (Hagenbach, 1822), which has mostly 
long- but also short-winged individuals, males are more mobile 
than females (Poniatowski and Fartmann 2011). While there may 
be other reasons for this difference, the bioacoustic environment is 
thought to be the most important factor (Berggren 2005). 
Although the trend of our SSD results indicated that the Zigana 
Mountain population may be more dimorphic than the southern 
populations, this difference was not statistically significant. Some in- 


sect and many other animal species have been reported to have higher 
SSD in areas with higher humidity and rainfall (Stillwell et al. 2007). 
However, wider sampling is needed to determine why SSD is weak 
in JA and to understand its relationships with the ecological factors 
affecting populations in different geographical areas. It is possible 
that no specific life-history traits emerged because the sampled indi- 
viduals of both sexes were in the same young age group in the labo- 
ratory and experienced similar feeding and growth environments. 


Bioacoustics and stridulatory file—As in most Orthoptera groups, 
the male calling song and stridulatory organ morphology of 
Isophya species have species-specific characteristics (Zhantiev and 
Dubrovin 1977, Heller 1988, Heller et al. 2004, Sevgili 2004, 
Zhantiev et al. 2017) that are particularly important in identify- 
ing cryptic species (Sevgili et al. 2006, 2012, Iorgu 2012, Sz6vényi 
et al. 2012, Sevgili 2018, 2020). Although the male calling song 
in Orthoptera has largely species-specific features, song temporal 
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Fig. 10. Distribution of carbon atom numbers in cuticular hydro- 
carbon profiles of males and females across all populations studied. 


structure, sound signal transmission, and hearing are also affected 
by environmental conditions, such as ambient temperature and 
vegetation structure, and anthropogenic factors, such as noise and 
habitat fragmentation (Hall and Robinson 2021). Isophya species 
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generally exhibit a fragmented habitat distribution pattern, and 
intraspecific variation is inevitably common given their extremely 
limited mobility. For example, I. modestior has a wide and frag- 
mented distribution in the Balkan Peninsula and has variations 
in the duration of song, the number of stridulatory organ teeth, 
and the impulse of the main syllable, while populations in the 
northern (Pannonian) and southern (Balkan) regions form two 
groups (Ivkovic¢ et al. 2022). Another study using two mitochon- 
drial (COI and 12S) and two nuclear (ITS2 and H3) molecular 
markers, the presence of two clades identified based on acoustic 
data was supported (Ivkovic¢ et al. 2023). However, it was shown 
that the Pannonian-North clade, which includes Austria, Serbia, 
Hungary, and Slovenia, also has subgroups. 

We found that both PC1 and PC2 are determined mainly by 
particular temporal features of the song and the structure of the 
stridulatory organ (Figs 7, 8). This is important in revealing vari- 
ations in secondary sexual characteristics in the partially isolated 
local populations of the same species. The observed environ- 
mental variations inferred from remote sensing data, including 
temperature, humidity, precipitation, NDVI, NDMI, and SWIR, 
among different localities suggest the exposure of local popula- 
tions to distinct selective pressures, thereby indicating potential 
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Fig. 11. Principal component analysis (PCA) plot of cuticular hydrocarbon (CHC) phenotypes of males and females from populations 
in different localities (A). Sexually dimorphic CHC differentiation in the species when all locations are included (B). The ellipses in the 
figure are the confidence ellipses around the group means in the principal components space of the function. 
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adaptation processes. In tettigonids, increased ambient tempera- 
ture significantly differentiates features such as calling song du- 
ration, repetition rate, and number of impulses in each syllables 
(Arias et al. 2012, Cooper et al. 2013, Cusano et al. 2016). These 
studies are important in showing that some environmental factors 
induce variation in some parameters of calling songs and point 
to possible causes of differences between at least some of the JA 
populations. For example, Cooper et al. (2013) reported regional 
differences for Tettigonia viridissima (Linnaeus, 1758) in both mor- 
phology and bioacoustic parameters (e.g., minimum stridulation 
frequency). However, in another study, T. viridissima populations 
were found to exhibit relatively uniform song patterns over a wide 
morphological range (Grzywacz et al. 2017). It is also worth not- 
ing that this species has long wings and is highly mobile. 


Cuticular hydrocarbon profiles. —Regarding variation among local 
populations, CHC profiles did not show significant variation in 
the morphological structures examined for both sexes, as well as 
in the stridulatory organ located of the left wing and bioacoustic 
characters of males. 

Several studies on CHC profiles have been conducted on vari- 
ous insect orders, such as Coleoptera (Niogret et al. 2019), For- 
micidae (Hymenoptera) (Martin et al. 2008, Morrison and Witte 
2011), and Diptera (Moore et al. 2021, Kula et al. 2023). The 
findings suggest that CHC profile differentiation can be used as 
a species-level chemotaxonomic tool. Although CHCs are highly 
stable (Drijfhout 2010), CHC chemical profiles vary significantly 
across geographical regions in some insect species, while in oth- 
ers, the differentiation is weak or absent (Suvanto et al. 2000, 
Moore et al. 2022). In addition, CHC profiles may vary to protect 
different populations of the same species from drying out or from 
cold/heat due to abiotic environmental variation (such as habi- 
tats having different climatic characteristics) (Veltsos et al. 2012, 
Karpati et al. 2023). CHC profiles may also vary with the age and 
seasonal morphs of the insect (Karpati et al. 2023). In our study, 
in which the JA individuals were of similar age, there were no 
significant differences in CHC profile compositions between all 
local populations except for one. 

Although CHCs can be useful in distinguishing between spe- 
cies, they may not differentiate as strongly as morphological, bio- 
acoustic, and genetic characters (Morrison and Witte 2011, Kather 
and Martin 2012, Veltsos et al. 2012, Dos Santos and do Nascimen- 
to 2015, Karpati et al. 2023). On the other hand, close relatives and 
sympatric species of the meadow grasshoppers Chorthippus (Glypto- 
bothrus) biguttulus (Linnaeus, 1758) and C. (Glyptobothrus) mollis 
(Charpentier, 1825) show clear differences in the position of the 
first methyl group in major methyl-branched CHCs (Finck et al. 
2016). Through an analysis of the molecular (COI gene) and CHC 
profiles of six Acridid (Orthoptera) species, Sofrane et al. (2022) 
demonstrated interspecific differences using both molecular and 
biochemical (weaker) data. Since pre-mating communication in 
Isophya species occurs through male call signals and short female 
responses, pheromonal communication, possibly associated with 
the CHC profile, may not develop strongly. In Isophya, male calling 
songs exhibit clear distinctions between species. 

Studies on Isophya and Tettigonia bushcrickets have reported 
the molecular data largely matching the bioacoustic and mor- 
phological characters (but see; Grzywacz-Gibala et al. 2010, Cho- 
banov et al. 2017, Grzywacz et al. 2017). In our study, there was 
a clear SSD in CHC composition in the JA individuals collected 
from each location. Similar findings have been reported in many 
insect species, such as the Australian bushcricket Kawanaphila 
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nartee Rentz, 1993. These CHC profile differences can be used to 
indicate relatedness, mate quality, body size, or fertility (Hare et 
al. 2022). For the first three components of the CHC profiles, we 
found significant differences between the Vauk-Zigana moun- 
tain populations and the others. Genes regulating melanization, 
which provides resistance to water loss and desiccation in insects, 
can also affect CHC profiles (Lamb et al. 2020). For example, it 
has been shown in Isophya rizeensis that melanization, which ex- 
hibits an altitude/temperature dependent change, is also affect- 
ed by solar radiation and vegetation height in the area (Kuyucu 
et al. 2018). Differences in NDVI between the studied regions do 
not seem to be related to differences in IA CHC profiles. Given 
the distinct environmental conditions on Mount Vauk, includ- 
ing temperature, relative humidity, and average precipitation, 
compared to other regions, these factors likely contribute to the 
notable divergence in CHC profiles observed in the Vauk popu- 
lation compared to other local populations, particularly regard- 
ing PC1 (Table 1). Kota et al. (2021) reported differences in both 
sexual dimorphism and the effect of geographic isolation on the 
CHC profiles of different island populations of the Pacific field 
cricket Teleogryllus (Teleogryllus) oceanicus (Le Guillou, 1841), al- 
though these differences were not significantly affected by local, 
environmental, and climatic variables. Similarly, Dupraz et al. 
(2022) found different CHC profiles for separate populations of 
the seabird tick Ixodes uriae White, 1852, and Kula et al. (2023) 
found Calliphora vicina (Robineau-Desvoidy, 1830) (Diptera) to 
have a distinct CHC profile in European (Germany and England) 
and Turkish populations. 

Finck et al. (2016) detected 34 different peaks with 25-33 
carbon atom components in the meadow grasshoppers Chorthip- 
pus biguttulus and C. mollis. On the other hand, greater diversity 
in chain length and carbon atom number (males: 15-38; female: 
22-38) was found in JA. Hare et al. (2022) identified 23 CHC 
components with chain lengths ranging from 14 to 39 in Kawa- 
naphila nartee (Orthoptera: Tettigoniidae). Studying the flesh fly 
Sarcophaga Meigen, 1826 from Diptera, Moore et al. (2021) report- 
ed a profile containing alkanes, methyl and di-methyl alkanes, and 
alkenes with chain lengths varying between C23-C33, which dif- 
fered between species and sexes within each species. Overall, these 
results indicate a wide variation in CHC profiles both within and 
between sexes in the genus Isophya. 


Conclusion 


Our findings reveal a high degree of uniformity across morpho- 
logical, bioacoustic, and CHC profiles in the studied populations of 
IA. However, noteworthy intraspecific variations were observed in 
certain localities, particularly regarding male stridulatory organ and 
bioacoustic parameters. Notably, the Vauk population exhibited 
unique CHC profiles compared to the other populations. This vari- 
ability is evident across different geographic regions, although it did 
not consistently follow the same pattern for all studied characteristics. 

Furthermore, our study provides preliminary insights into in- 
traspecific sexual dimorphism and sexual selection in JA. However, 
further research is necessary to fully understand these variations 
across different IA populations. 

Thus, this study contributes to a deeper understanding of intraspe- 
cific variability within IA. As the habitats of narrow-range endemic 
species such as IA face threats from human impacts, such as fragmen- 
tation, degradation, and even elimination, understanding the varia- 
tions between local populations becomes crucial for the conservation 
of these ecologically diverse endemic populations in Anatolia. 
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Supplementary material 1 


Author: Ebru Kiran Ozdemir, Hasan Sevgili, Emine Bagdatli 

Data type: pdf 

Explanation note: fig. $1. Two males and a female with fresh sper- 
matophore of Isophya autumnalis and views of the species’ 
preferred habitats from the Poske and Zigana mountains. fig. 
S2. Body parts used for morphological measurements. fig. S3. 
SEM image of the stridulatory file of a male and morphometric 
measurements. fig. S4. Images plotting of normalized differenc- 
es vegetation index (NDVI), normalized differences moisture 
index (NDMI), and short-wave infrared (SWIR) for the local- 
ity from Zigana mountain. fig. $5. Images plotting of normal- 
ized differences vegetation index (NDVI), normalized differ- 
ences moisture index (NDMI), and short-wave infrared (SWIR) 
for the locality from Krom Valley. fig. S6. Images plotting of 
normalized differences vegetation index (NDVI), normalized 
differences moisture index (NDMI), and short-wave infrared 
(SWIR) for the locality from Vauk mountain. fig. S7. Images 
plotting of normalized differences vegetation index (NDVI), 
normalized differences moisture index (NDMI), and short- 
wave infrared (SWIR) for the locality from Pekiin mountain. fig. 
S8. Images plotting of normalized differences vegetation index 
(NDVI), normalized differences moisture index (NDMI), and 
short-wave infrared (SWIR) for the locality from Péske moun- 
tain. fig. S9. Comparison of average normalized differences veg- 
etation index (NDVI) data for July and August according to lo- 
calities. fig. S10. Comparison of average normalized differences 
moisture index (NDMI) data for July and August according to 
localities. fig. $11. Comparison of average short-wave infrared 
(SWIR) data for July and August according to localities. 

Copyright notice: This dataset is made available under the Open 
Database License __(http://opendatacommons.org/licenses/ 
odbl/1.0/). The Open Database License (ODbL) is a license 
agreement intended to allow users to freely share, modify, and 
use this Dataset while maintaining this same freedom for others, 
provided that the original source and author(s) are credited. 

Link: https://doi.org/10.3897/jor.33.118937.suppl1 
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Supplementary material 2 


Author: Ebru Kiran Ozdemir, Hasan Sevgili, Emine Bagdatl 

Data type: docx 

Explanation note: table $1. Monthly and annual temperature and 
relative humidity data for the studied region and its vicinity. 
No data available for Pekun mountain. table S2. Results of 
Tukey HSD statistics of some remote sending data between lo- 
calities. NDVI: Normalized difference vegetation index. NDMI: 
Normalized differences moisture index. SWIR: short-wave in- 
frared reflectance. table S3. Measurements of morphology and 
some bioacoustic parameters in I. autumnalis. 

Copyright notice: This dataset is made available under the Open 
Database License _(http://opendatacommons.org/licenses/ 
odbl/1.0/). The Open Database License (ODbL) is a license 
agreement intended to allow users to freely share, modify, and 
use this Dataset while maintaining this same freedom for others, 
provided that the original source and author(s) are credited. 

Link: https://doi.org/10.3897/jor.33.118937.suppl2 
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